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Quantum Fluctuations of Chirality in One-Dimensional Spin-| Multiferroics: 
Gapless Dielectric Response from Phasons and Chiral Solitons 
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We present a quantum theory for one-dimensional spin- 1/2 multiferroics, where the vec- 
tor spin chirality couples with the electric polarization. Based on exact diagonalization and 
bosonization, it is shown that quantum fluctuations appreciably reduce the chiral ordering am- 
plitude and the associated ferroelectric polarization. This yields nearly coUinear spin correla- 
tions in short-range scales, in qualitative agreement with recent neutron scattering experiments. 
There appear gapless chirality excitations described by phasons and new solitons, which can 
be experimentally verified from the low-energy dielectric response. 
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When the chiral symmetry is spontaneously bro- 
ken by magnetic interactions in Mott insulators, a ferro- 
electric polarization appears through the spin-orbit cou- 
pling. This new prototype of multiferroic behavior has 
been discovered in a spin-2 helimagnet ThMnOa,^-* and 
attracted a current great interest for both its funda- 
mental importance and its potential application to an 
electrical control of spins. ^■'^^ This magnetoelectric cou- 
pling between the vector spin chirality and the polariza- 
tion'^"^) is ubiquitous in Mott insulators. ''^ Of our par- 
ticular interest is recently discovered multiferroic behav- 
ior in one-dimensional (ID) frustrated spin-1/2 magnets, 
LiCuV04^"") and LiCuzOs.^^""' These have opened an 
intriguing issue of strong quantum fluctuations in multi- 
ferroics due to the low dimensionality and the spin-1/2 
nature, as suggested by recent neutron scattering exper- 
iments."' 

Another important aspect in multiferroics is that 
low-frequency magnetic excitations can be probed 
from the dielectric functions e"(ijj) through the mag- 
netoelectric coupling, ^^"^^^ as actually measured for 
i?Mn03.^^'^°) Namely, a local flip of the chirality in- 
duces the charge dynamics. Particularly, in the ID spin- 
1/2 multiferroics, a pair of such local defects may prop- 
agate as new elementary excitations. This is reminiscent 
of charged solitons in ID band insulators, e.g., polyacety- 
lene.^"'^-' However, the chirality-induced charge dynamics 
in the ID multiferroics has a much smaller energy scale 
and may produce the low-frequency dielectric response. 

In this Letter, we for the first time reveal the novel 
quantum nature of ID spin-1/2 multiferroics, by ex- 
amining a simple but realistic frustrated spin model. 
It is shown that the chiral phase can appear near the 
SU (2)-symmetric case, which explains the experimen- 
tally observed ferroelectricity. Significantly, this multifer- 
roic state is characterized by (i) a tiny amplitude of the 



the chiral long-range order (LRO), (ii) almost coUinear 
spin correlations, and (iii) gapless chiral solitons in the 
low-frequency dielectric response. 

The minimal model for the ID spin-1/2 multifer- 
roics, especially for LiCuV04,^^) is given by a simple spin 
Hamiltonian with nearest- and second-nearest-neighbour 
exchange couplings, Ji and J2: 
2 

^ = E E"^" ■ ^^■+" + ^ i)5|^i+„] • (1) 

n=l j 

Here, 5'^ represents the spin operator on the site 




Fig. 1. (Color online) The intensity plot of the spin Drude weight 
Ds/J2 for the model (1) with L = 28. The black broken line 
at A 1 shows a chiral phase boundary obtained in the 
scaling analysis. The red lines around IJ1I/J2 ~ 3 - 4 are the 
phase boundaries between TLLs and dimer phases determined in 
refs. 28 and 29. At Ji = 0, there is no chiral order. A dark region 
with Ji < contains at least two phases; for small A, a dimer 
order appears^^^ with a unit lH) -|- |J,t) on the nearest-neighbor 
bonds, while for A Ri 1, another gapped phase is expected to 
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j (= 1,...,L) and A(< 1) is an easy-plane exchange 
anisotropy. In the classical limit, a chiral ground state is 

realized when | J1I/J2 < 4, having a nonzero expectation 
value = {O^} of the vector spin chirality 



x,y,z. 



(2) 



In the quantum case with spin-1/2, the existence of 

a chiral phase with an incommensurate gaplcss excita- 
tion has been proposed based on the bosonization com- 
bined with a mean field theory,^^' and has been con- 
firmed numerically aroimd the XY case A « 0.^"'' On 
the other hand, in the 5'C/(2)-symmetric case A = 1, 
the chiral ordering is prohibited, and valence bond solids 
(VBS) are stabilized. The intermediate regime of 
the anisotropy < A < 1 has not been clarified yet. We 
revisit this issue by calculating the spin Drude weight Ds 
by Lanczos diagonalization,^^) to see whether the spin 
excitation is gapless or not; see Fig. 1. We find a re- 
gion with moderate Ds/J2 for IJ1I/J2 < 1, which can 
be assigned to the gaplcss chiral phase found in rcfs. 
23 and 24. When increasing IJ1I/J2 from zero in the 
XY case A = 0, Ds/J2 starts to decrease rapidly at 
\Ji\/J2 ~ 0.9, in reasonable agreement with the previ- 
ous numerical results on the phase boundary between 
the chiral and dimer phases. When increasing A 
from zero for IJ1I/J2 ^ 1, Ds/J2 decreases only gradu- 
ally, which suggests that the gapless chiral phase might 
extend to the vicinity of A = 1. On the other hand, it is 
known^^) that the gap associated with the VBS states for 
A = 1 and | J1I/J2 ^ 1 is too small to detect in numerical 
analyses of small clusters. 

To correctly determine the phase diagram for small 
IJ1I/J2, we present a scaling analysis based on the 
bosonization. Wc start from the bosonized effective 
Hamiltonian of cq. (1) formulated in ref. 23 for | J1I/J2 ^ 
1, which consists of two Tomonaga-Luttingcr liquid 
(TLL) parts 7i±[<i>±,6±] and interaction terms TYint- 
Here ("I>y, 8,^) is the pair of dual boson fields. The TLL 
parameters of the sectors H± are obtained as 



K4 



K 



2-KV \\J2 




(3) 



where K and v are the TLL parameter and the spin- 
wave velocity for Ji = 0, respectively. With the easy- 
plane anisotropy A < 1, two interaction terms Oc = 
(90+) sin(\/47r9-) with the scaling dimension Ac = 
1 -I- 1/K- and Od = cos(A/47r$+) cos(V47r0-) with the 
dimension A^ = K+ + 1/K- can be relevant for Ji > 0, 
while only Oc can for Ji < 0. Here, a relevant C'c(rf) can 
yield a gapless chiral (singlet-dimcr) order. Therefore, 
we obtain the following criterion: for Ji > 0, the phase 
boundary between dimer and chiral phases is fixed by 
the relation Ac = A^, while for Ji < 0, the chiral or- 
der appears under the condition Ac < 2. Then the chiral 



phase boundary is given by 

A « 1 - (Ji/ J2)V(2^') + • • • (4) 

for IJ1I/J2 1, as indicated by the black broken line 
in Fig. 1. Namely, a small easy-plane anisotropy, which 
is realistic for spin-1/2 systems, induces the chiral order. 
We also calculate the one- loop RG corrections to K±, 
which suggest that the chiral phase is robust for Ji < 
while it is reduced by an expansion of the dimer phase for 
Ji > 0. This tendency is in accord with the slower decay 
of Ds/J2 on the ferromagnetic side Ji < in Fig. 1. 
Henceforth, we focus on the realistic case A w 1 for the 
cuprates. 

Now let us investigate properties of the chiral ground 
state, i.e., quantum multiferroics. The imiform chiral- 
ity K,^ produces the electric polarization through the 
inverse Dzyaloshinskii-Moriya (DM) interaction,^"'^' i.e., 

cx ect/3xCf , where eap-^ is the fully antisymmetric ten- 
sor. Then it also generates the lattice modulation due to 
the electrostatic force. We take into account this coupling 
between spins and transverse optical phonons through 
the DM interaction."^^' Integrating out phonons gives a 
biquadratic DM interaction. 



-ffBDM = -V^'^{Sj X Sj+iY{Sj+2 X Sj+i) 



(5) 



Here, longer-range terms have been neglected for simplic- 
ity. This term yields a gapless chiral phase for quantum 
spin ladders"^^' and a spin cholesterics for classical spin 
systems.^^'^^' 

First, the chiral order parameter calculated from 
= ^/ {{^kY) is shown in Fig. 2. The finite-size re- 
sults are extrapolated to L ^ 00 by the alternating e- 
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Fig. 2. The left panel shows the chiral order parameter = 
-JI-SPlY) for J1/J2 = -0.1 and A = 0.9 as a function of 
\''zl J'2.- Extrapolation is done from the sequence of data for L = 
16, 20, 24, 28 and 32, using the alternating e-algorithm.^*) On the 
right panel, the extrapolated values (placed at "1/L = 0") are 
plotted together with the finite-size data. Broken lines are guides 
to eyes. 
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algorithm.'^'*' We take a small | J1I/J2 to assure that the 
system is inside the chiral phase even at = and 
to make the incommensurate wave vector Q sufficiently 
close to 7r/2. Remarkably, when T4/J2 < 2 for A = 0.9, 
is significantly suppressed from the classical value 
S*^ = 1/4, indicating appreciable quantum fluctuations of 
the chirality near the S'?7(2)-symmetric case. This is con- 
sistent with the bosonization analysis combined with the 
mean-field theory, which gives ~ \Ji/ J2\^^'^^^~^^ 
with K- — 1 ~ ^/\ — A near A = 1. On the other hand, 
with increasing T4, k,' starts to increase rapidly around 
Vz/ J2 ~ 2 and approaches a constant I/tt, which is the 
maximal value for the spin-1/2 case. 

The weak chirality for 14 /J2 J; 2 has a remark- 
able consequence on the spin correlations. Figure 3 (a) 
shows the equal-time in-plane spin correlation, S'^ {q) — 
(S+SZ^) with = Ej5"e-*«V\/Z for V4/J2 = 2. To 
break the inversion symmetry in finite-size systems and 
to have a chiral order parameter comparable to the value 
extrapolated to L 00 in Fig. 2, a flux = — O.Stt is 
inserted through the periodic chain as a mean field. Fig- 
ure 3 (a) shows two comparable peaks at q = ±Q and a 
tiny asymmetry under the inversion q —q, reflecting 
a small chirality. According to the bosonization analysis 
of the gapless chiral phase, ^'^^ only the peak at q — Q 
eventually diverges in the thermodynamic limit, due to 
a spiral quasi-LRO. The other peak at g = —Q observed 
for Vzl J2 = 2 converges to a large but finite value. Nev- 
ertheless, except in the close vicinity ai q — iQ, the 
asymmetry is small for Vz/ J2 ^ 2. To understand the 
physical meaning of these results, it is useful to intro- 
duce 



a{q) 



^/S+^)-^S+-{-q) 



(6) 



^fS+^+^S+-{-q) 
for each q component, which represents the ratio of the 




Fig. 3. (Color online) (a) The spin structure factor cS+ (g) for 
J1/J2 = -0.1, A = 0.9, and Vz/J2 = 2 with a flux = -O.Stt. 
The insets show the </)-dependence of /t^ = {0%). (b) The elliptic 
ratio a{q) for each q component of the spiral correlations, for the 
same J1/J2 and A and for Vz/J2 = 2 and 4. 



minor to major axes of ellipse. When a = (a = 1), the 
spin correlation is collinear (a circular spiral). The diver- 
gence of (Q) indicates a{Q) = 1 in the thermody- 
namic limit. On the other hand, away from q = Q ^ 7r/2, 
the elliptic ratio a{q) converges to a rather small value, 
particularly for smaller V2/J2, as shown in Fig. 3 (b). 
Therefore, the spin correlations are elliptic, even nearly 
collinear, at short-range scales and only the quasi-LRO 
component is characterized by a circular spiral with a 
pitch Q. The circular to collinear crossover in the length 
scale also appears in the time scale. Namely, the dynam- 
ical structure factor diverges algebraically as 
'-^ and our bosonization analysis shows 
that an energy gap Eg ^ J2\Ji/ J2\^^ opens 
at (7 = — Q, which is tiny near the S'C/(2)-symmetric 
case. Thus, 5^^(g,a;) shows a circular spiral only for 
uj, v\\q\ — Ql ^ Sg, and is nearly collinear otherwise. 

This provides a key to understanding puzzling results 
from a recent elastic polarized-neutron scattering exper- 
iment in LiCu202:^''-' the elastic elliptic ratio sa 9 - 20% 
indicating an elliptic spiral, and S^^{Q,Q) w Syy{Q,0) 
indicating a circular spiral, if one interprets that the ob- 
served intensity originates solely from the genuine Bragg 
peak due to a 3D magnetic LRO. This problem can be re- 
solved, at least qualitatively, if the fine but flnitc energy- 
momentum resolution allows an additional contribution 
from low-energy spin fluctuations to the elastic intensity 
(see also a related argument in ref. 18). Actually, accord- 
ing to our analysis of the weak-chiral region Vzl J2 ^ 2, 
there exist large spin fluctuations outside the narrow win- 
dow w, ujlgl — Ql ^ £g, which show small a but retain 
gxx _ gyy gygjj high uj. Detailed inelastic polarized 
neutron scattering experiments are required to test this 
possibility and to reveal nontrivial effects of a weak 3D 
coupling on the spin dynamics. 

Next, we discuss dynamical in-plane and out-of- 
plane chiral correlations, K'' and ^ respectively; 



K"'{q,uj) 



- 5(GS|0;?(q)- 

TT UJ - 



1 



n + iij 



(7) 

where 77 ^ -f and 0'^{q) = L-^ 'Zj e""'^ [Sj x S'j+i)" is 
the Fourier component of the chirality. There appear the 
following contributions to dielectric functions e°'^{q,uj) 
through the magnetoelectric coupling: !3e^^(g,a;) cx 
K^{q,Lu) and 3e^^(g, w) cx Ky{q,uj). We have performed 
numerical continued-fraction analyses of Lanczos results 
on K°'(q,uj) at q = 0,7r. Then, it is found that Ky{G,uj), 
which was argued in the context of electromagnons,^^-* 
and Ky{Tr,uj) are gapful. On the contrary, both K^{0,uj) 
and K^{tt,uj) are gapless, as we discuss below. To be 
explicit, we show the numerical results on if'^(0,a;) and 
K^{tt, ui) in Figs. 4 (a) and (b), respectively, for the same 
J1/J2 and A as used in Fig. 2 but for Vz/J2 = 3. In the 
insets, we clarify finite-size effects on the energy levels 
of the first excited states (A) and a higher-energy state 
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(B) with large intensity. They indicate that both chiral- 
ity excitations K^{0, iv) and K^{-k, lu) eventually become 
gapless in the thermodynamic limit, leading to a gap- 
less dielectric response even in a Mott insulator. This 
indicates an illuminating feature that this ID multifer- 
roic state possesses an antiferro-chiral (and thus antifer- 
roelectric) quasi-LRO, in addition to the uniform chiral 
LRO. 

Now we clarify the origin of the above gapless chi- 
rality dynamics using the bosonized expressions of the 
chirality operators in the gapless chiral phase: 

— [ci sin(V4^e_) + C2ag«;^(a,e+)2 



/dx 
- — csk"" cos(V47r$^ 
Lao 



(8a) 
(8b) 



where Ci_3 are 0(1) constants and is a lattice spac- 
ing. Note that excitations associated with (<i>_,G_) are 
gapful because of the chiral LRO, while those associated 
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Fig. 4. (Color online) Dynamical in-plane chirality correlation 
function K^{q,w) at (a) q = and (b) q = it, respectively, for 
J1/J2 = —0.1, A = 0.9, and Vz/J2 = 3. The shaded areas show 
the spectra for L = 28 drawn with r} = 0.2 J2 (see eq. (7)). In (a), 
the pseudo ground state is removed since it does not contribute 
to the dynamics in the thermodynamic limit. The symbols show 
the matrix elements LM^ ^ = L''\(\\0^{q)^GS)\^ at oj = ujx, 
where |A) denotes an eigenstate with energy ljx measured from 
the (true) ground state. For L = 28, they are shown for all A's; for 
smaller L, they are shown only for the low-energy levels depicted 
as A and B. The insets show the extrapolations of the A and B 
levels into L — » 00. (c-e) Schematic pictures on the creation and 
propagation of the chiral soliton pair in the background of quasi- 
long-range spiral order. The symbol L (R) denotes a positive or 
left-handed (negative or right-handed) chirality at each bond. 



with (<&+, 0+) are gapless since the spin rotational sym- 
metry is conserved. Let us begin with the uniform part. 
The first term in eq. (8a) gives rise to gapful spinon-pair 
excitation in K^{0,u!), which is assigned to the shaded 
spectrum above the energy level B in Fig. 4 (a). On the 
other hand, the second term in eq. (8a) gives a nonsingu- 
lar gapless contributions to K^{0, w), which explains the 
spectrum from A in Fig. 4 (a). This gapless mode turns 
to two-phason excitations if a 3D incommensurate spiral 
magnetic LRO appears due to a weak inter-chain cou- 
pling. In realistic systems containing defects that act as 
pinning centers, single phason excitation can also appear 
in ii:^(0,w).35) 

Similarly, the term cos(-\/47r*I'+) in the staggered 
part Of.{TT) also yields a gapless mode, which leads to 
power-law behavior K^{tt,uj) ~ see Fig 4 (b). 

The physical origin is understood by noticing the follow- 



ing properties: (i) cos(v47r$+) translates 8+ in Sj oc 
gi\/7Fe+ |-|y g^jjj j-j^-j Jdx cos(\/47r'i'+) coincides with 
the leading term in the bosonized form of J2j SjSj_^_i, 
which rotates the two neighbouring spins by tt about the 
z axis. Accordingly, the gapless excitations in K''{Tr,uj) 
consist of two solitons, as schematically shown in Fig. 4 
(c)-(e): An operation of cos(\/47r<l)+) transforms a state 
from (c) to (d). It creates a pair of kinks, which con- 
tain, for instance, a negative (i?) chirality in the posi- 
tive (L) chiral background. Hence, they induce the op- 
posite transverse displacement of electric charge in com- 
parison with the ferroelectric background. Then, via the 
exchange of two neighboring chiral spin pairs denoted by 
circles, the kinks propagate without changing the sign of 
the chirality of each domain, as shown in (e). This soliton 
results from a chiral pairing of two spinous, in contrast to 
the gapless charge-neutral spinon in ID spin- 1/2 antifer- 
romagnets. The transverse nature of the charge displace- 
ment also contrasts with the longitudinal nature of the 
well-known charge soliton in the polyacetylene.^^^ In the 
ID multiferroics like LiCuV04 where the unit cell is dou- 
bled, the chiral solitons observed in the staggered part 
can also appear in the uniform part. If the U{1) symme- 
try is weakly broken by spin anisotropy or an in-plane 
magnetic field, the soliton can acquire a small gap. It will 
be interesting to look for this soliton-pair excitation by 
low-frequency electrical or optical probes. 

In real materials, an inter-chain coupling J' causes 
the 3D spiral LRO. However, a low transition tempera- 
ture Tn « 2K in LiCuV04^""^ suggests that the effect 
of J' is small and should not alter our main conclusions, 
particularly on the quantum dynamics. 
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